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Inferring the dynamics of stellar streams via distance gradients 
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ABSTRACT 

We present a simple result in which the distance gradient along a stream can be used to derive 
the transverse velocity (i.e. proper motion) along it, if the line-of-sight velocity is also known. 
We show its application to a mock orbit to illustrate its validity and usage. For less extended 
objects, such as globular clusters and satellite galaxies being tidally disrupted, the same result 
can be applied in its small-angle approximation. The procedure does not rely on energy or 
angular momentum conservation and hence does not require a Galactic model in order to 
deduce the local velocity vector of the stream. 
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1 INTRODUCTION 

In the last few years, our view of the Milky Way's stellar halo has 
been altered by the pioneering work of the Sloan Digital Sky Sur- 
vey (SDSS) and continued analyses of its data; it is now abundantly 
clear that the outskirts of the Galaxy provide the prime locations 
in which to search for signs of ongoing accretion events onto our 
Galaxy, as well as the remnants of such events, more progressed 
toward their final dissolution into the Milky Way teelokurov et al.l 
l200fl ). Current survey data that enable such studies are mainly pho- 
tometric, while spectroscopic follow-up is often necessary to check, 
via its kinematics, whether a photometrically selected stellar group- 
, ing is truly an association. Many stellar streams originating from 
■ either globular clusters or dwarf galaxies have been identified in 
the Galactic halo through such investigations. The form er can be 
,. exemplified by the globular cluster Palomar 5 ( Odenkir chen et al.l 
' l200lll2003h . while the most promin ent example of the latter is the 
Sagittarius stream jlbata et allll994[) . 

Much of the recent attempts to understand the dynamics of 
these objects, and thus also their ultimate fate, have come from 
their modelling via numerical simulations that can now foster mil- 
lions of particles in a study of a single stellar system. These works 
illustrate how the tidal distortion and subsequently resulting disrup- 
tion lead to orbits of constituent stars that dev iate from that of the 
bound pare ntal struct ure (e.g.llb ata et al. 1997; Dehnen et al. 2004; 
IJohnston et al, 2005i ; lchoi et al..2007 ; Montuori et al. 20 07). Fully 
live numerical studies also allow for the possibility of studying dis- 
rupting systems under the auspices of dynamical friction. Simu- 
lations are therefore key to furthering our understanding of these 
dynamically evolving stellar substructures, with improvements and 
advancement continually being made on the incorporation of neces- 
sary physics in the correct ratios to model better the systems under 
study. 
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It is, however, also interesting to reflect on the extent to which 
simple geometric or analytic arguments can aid in these quests. Pre- 
viously, we showed that the radial velocity gradient along a stream 
can be used to derive th e transverse velocity (i.e . proper motion) 
as a f unction of distance jjin & Lvnden-Belll2007l ; see also lsinnevl 
l2008l) . which in turn allows for the determination of orbits for these 
streams. The success of this method relies on the radial velocity 
gradient being well determined, and Binney showed its successful 
application to Palomar 5, by recovering its observationally known 
properties. The procedure presented in this letter is a related and 
complementary one, applicable when the distance information for 
stars in a stream is more abundant (or better constrained) than that 
for radial velocities. In the following sections, we first present a 
simple result that can be used to calculate the transverse veloc- 
ity for a given point along a stellar stream, using knowledge of 
the distance gradient and the radial velocity there. An application 
to a mock orbit is then used to illustrate its usage. We also pro- 
vide a simple geometric visualisation of how the same result can 
be reached for less extended objects, such as dwarf galaxies that 
are being tidally disrupted. 



2 CONSTRAINING TRANSVERSE MOTION THROUGH 
DISTANCE GRADIENTS 

2.1 The method 

Let dl be the heliocentric position vector to a particular point on the 
stream, with unit line-of-sight vector 1 and heliocentric distance d, 
so that its line-of-sight velocity, relative to the Galactic Standard of 
Rest, is given by vi — v.l. Differentiating the position vector with 
respect to time, t, and using the chain rule results in the following: 

dd- ,di 
V = — 1 + d— 

dt dt 



^^1 + d—^ 

dx dt dx dt 
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Vs dd- 
d dx 



(1) 



where x is the angle along the stream, measured from some fiducial 
point, and s is the unit vector along the apparent direction of mo- 
tion of the stream at 1, such that s — dl/dx and dx/dt — v.s/d = 
Vs/d from the usual equation for the transverse velocity. One pos- 
sible m ethod for calculating § is provided by I jin & Lvnden-Belll 
( l2008h . By noting that the full velocity vector at any point along a 
stream is given by v = vil + VsS, we obtain the following relation: 



dvi 



(2) 



This result can also be reached by taking the derivative of the po- 
sition vector with respect to x^ rather than t. We note that equa- 
tion (5) of Binney (2008) gives the same result, by identifying the 
time derivative of the angle as v^/d. 

While equation ^ does not require the value of the Galactic 
potential at the location of interest as an input, it assumes that the 
stars of the stream, whose distances are used as an ingredient of 
the method, are part of the same moving group and that the stellar 
trajectories follow the observed stream. Under this assumption, the 
method provides a means of calculating the stream's local trans- 
verse velocity, thereby allowing the determination of all compo- 
nents of the velocity vector. In an ideal case, one would have both 
distance and kinematic data on the stars. However, if both distances 
and line-of-sight velocities are known, then we could just as easily 
calculate the transverse velocity through a method that relies on the 
gradient in vi ( Jin & Lynden-Bell 2007). The result presented in 
this letter is therefore most useful when the distance gradient along 
a stream is well determined, and for relatively small uncertainties 
in the stellar distances (e.g. for RRLyrae or BHB stars), but when 
one's knowledge of the kinematics is more limited. On the other 
hand, if no radial velocity information is available for the stars, one 
can still determine the ratio of the two velocity components. Hence 
the family of orbits to be derived for the stream need only be a 
function of vi . 

Given the simple form of equation the error in the derived 
transverse velocity is calculated very straight-forwardly through the 
quadrature sum of the fractional errors in the other quantities: 



+ 



Ad 
~d 



1/2 



(3) 



where Q — dd/dx and Ax gives the magnitude of the error, or 
uncertainty, in quantity x. This warns us that if any of the quanti- 
ties are very small, then its fractional error will significantly boost 
the error for Vs and that, in particular, a small distance gradient 
will cause the most damage, given its place in the denominator of 
equation l|2j- 



2.2 Example application 

Figure [T] shows the distance evolution along a mock orbit (in 
grey), generated using the Milky Way model adopted bv iPaczvnskil 
The exercise here is to use the distance gradient and a 
single radial velocity to deduce the transverse velocity at a given 
point on the stream, and then attempt to recover the input orbit. 
Within this orbit, we place ourselves at a point mid-way along 
the segment shown in Figure [T] at x = 8.0°. The distance gradi- 
ent is calculated by simply fitting a second-order polynomial (also 
shown in the same figure by the dot-dashed black line) and taking 
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Figure 1. Distance as a function of angle along a small segment of a mock 
stream (solid grey line), overplotted with a fit deduced using a second order 
polynomial, d = a + bx + cx^ (dot-dashed black line); the fit shown has 
a = 30.1, b = 0.509, c = 2.86 x 10"^. x = 8.0° denotes the location 
where we choose to calculate the distance gradient in order to deduce the 
transverse velocity, and hence the location for which the initial conditions 
for the orbit calculation are subsequently defined. 




Figure 2. Testing the recovery of the orbital paraineters: Galactic longitude 
(top panel), heliocentric distance (middle panel) and line-of-sight velocities 
relative to the Galactic Standard of Rest (bottom panel) are shown as func- 
tions of Galactic latitude. Original data froin the mock orbit are shown in 
dark grey (solid), whilst those for the recovered orbit are plotted in black 
(dashed). There are no significant differences between the two orbits over 
70°. Also shown in light grey are 147 orbits for which artificial uncertain- 
ties were added to the orbital starting point for the distance, line-of-sight 
velocity and distance gradient; values for each parameter are drawn from a 
Gaussian distribution centred on the value characterising the correctly re- 
covered orbit, with its Icr deviation set to 7% that of the distribution mean. 
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its derivative at tlie location of interesfl The heliocentric distance 
and line-of-sight velocity are 34.5kpc and 137.9 kms~^, respec- 
tively, and through the application of equation l|2}, we determine the 
transverse velocity to be 149.4kms~^. The directional vector § is 
found to be (-0.014, -0.254, -0.967) in Cartesian Gala ctic co- 
ordinates, using equation (16) of |jin & Lvnden-BeUl ^20081) on two 
locations, separated from the central point by 5°. These allow us 
to calculate the velocity vector, from which we then have the nec- 
essary 6-dimensional phase-space information to derive the orbit. 
The original orbit is properly recovered in this example, with the 
Cartesian velocity vector being (123.5, —93.4, —131.8) kms~^ in 
the deduced orbit at the position specified above, where it was 
(123.7, -93.0, -130.8) kms"^ in the original orbit. The small de- 
viations in the velocity vector components do not lead to any no- 
ticeable difference in any of the orbital parameters over 70°, as 
shown by the similarity between the solid grey and dashed black 
lines in Figure [2] Note that x ~ 8.0° in Figure [T] corresponds to 
a Galactic longitude and latitude of {I, b) — (—23.8, 5.3)° in Fig- 
ure (2] 

The panels of Figure[2]also show the results of a simple exer- 
cise, whereby the parameters d, vi and dd/dx at the orbital starting 
point were allowed to vary, in order to mimic and ascertain the pos- 
sible effects of observational uncertainties for a real stellar stream. 
Each of the three parameters were drawn from a Gaussian distribu- 
tion with the mean given by the initial conditions of the correctly 
recovered orbit, and whose la deviations were set to 7% of the 
value of the distribution mean. In practice, observational uncertain- 
ties are likely to be fractional for distances (~ 7% for RRLyrae 
and BHB stars) and absolute for velocities (~ 10 kms^^). For the 
parameters of our mock orbit, 7% of the mean value corresponds 
to a realistic uncertainty in each of the variables. The resulting set 
of 147 light grey lines in each of the panels therefore indicate the 
degree to which the recovered orbit might differ from the origi- 
nal in cases where we might expect the uncertainties for d, vi and 
dd/dx to sach be of the order of a few percent. As expected from 
this form of artificial 'degradation', the agreement of the sky posi- 
tions is remarkably good within ±20° of the orbital starting point, 
while the distance tends to suffer a constant offset and the recovery 
of the line-of-sight velocity is not as good. Unsurprisingly, devia- 
tions from the original orbit in all orbital parameters become larger 
with increasing separation from the orbital starting point, with devi- 
ations also manifesting themselves strongly near the turning point 
of the orbit. 

It is also necessary to add the usual cautionary note that 
streams do not follow exact orbits. However, the notion of quan- 
tifying their transverse velocity in this way provides a useful tool 
with which the motions of streams can be studied with limited kine- 
matic information. Figure [2] also highlights that, in general, one 
should preferentially apply the technique — observations allowing 
— every ^ 20° along the stream, so as to follow better the stream's 
properties and hence determine a realistic orbit. 




Figure 3. Schematic diagram to illustrate how equation can be reached 
in the .small-angle regime for objects with a more limited angular extent 
than a generic stellar stream. A denotes the position of the Sun, d is the 
heliocentric distance to the object at B, and EC gives the true extent of 
the object with the object's elongation being along the grey arrow. The 
difference in distance between the extreme ends of the object is given by 
Ad, where Ad << d. Equation (5) is recovered straight-forwardly when 
X < < 1 < < and the observed elongation of an object is assumed to be 
a simple projection of its true elongation onto the plane of the sky, with the 
orbit following the direction of elongation, such that tan ij> = fi/vs . 



galaxy becoming unbound and is thus no longer self-gravitating. 
Such a system would exhibit a large elongation and a radial ve- 
locity gradient along its extent. In these cases, one assumes that 
the observed elongation can be de-projected to give the true or- 
bital direction, whose angular difference to the sky projection, is 
given by the relation tan (f> = vi /vs. This implies that if the trans- 
verse velocity can be estimated through other means, then the vari- 
ation in distance acr oss the face of the galaxy can be easily deduced 
jMartin & Jii]l2009l) . Figure [3] illustrates how equation ^ may be 
reached in such cases, with the pictorial representation showing 
quantities in the plane of the orbit. Note that this figure is not to 
scale as, in reality, Ad/d « 1. In cases where x is small (or, 
more specifically, when x << 1 << 0), equation ^ is recovered 
by applying the small-angle approximation directly to Figure[3]and 
assuming that the true elongation (in the direction of the grey ar- 
row) follows the velocity vector. We also outline below the steps 
taken to reach the same answer when leaving the application of the 
approximation to the end; in this case, one can initially employ the 
sine and cosine rules to obtain exact relations between the various 
lengths and angles involved. 

Let L be the true extent of the object, denoted by BC in Fig- 
ure[3] Then: 



d 



L 



cos((;/> + x) sinx ' 

{d + Adf = d^ + L'^ + 2dL sin cf, 



(4) 



(5) 



2.3 Small-angle approximation 

The same relation for Vs can be reached in a straight-forward man- 
ner in the small-angle regime through simple geometric arguments, 
when the object of interest is, for instance, a disrupting dwarf 

^ The choice of location is reasonably random, apart from avoiding regions 
with little distance variation; this would lead to unreliable values for Vs as 
explained in the main text. 



and 

d"^ ^ {d + Adf + L"^ - 2{d + Ad)Lsm{(l) + x) ■ (6) 

By taking the difference between equations l[5) and ((SJ, and then 
using l|4j to eliminate L, one arrives at equation l|2j in the limit 
that X << 1 << and Ad « d. The latter assumption is only 
necessary in so far as the fact that the assumed linear elongation, L, 
should be a good representation of the direction of orbital motion 
at that location. 
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3 SUMMARY 

This letter presents a simple result that enables the calculation of 
the transverse velocity along a stream by using the gradient in dis- 
tance along it. If at least one radial velocity measurement is avail- 
able and the distance gradient there is known, the transverse ve- 
locity at this location can be calculated without the knowledge of 
the Galactic potential. In these cases, one then has all of the 6- 
dimensional phase-space information necessary to derive the orbit 
for the stream. For cases where the radial velocity is not known, one 
still gains knowledge of the ratio between the radial and transverse 
velocity components. It is therefore still possible to deduce fami- 
lies of orbits that are functions of the initial radial velocity alone. 
This result is particularly useful for the study of Galactic stellar 
streams, in cases where there is a lack of spectroscopic information, 
but when ample photometric data are available for stellar popula- 
tions with reliable distance measurements and with small associ- 
ated errors. With current and future large- sky photometric surveys 
— such as the SDSS and Pan-STARRS — having the potential to 
find a multitude of stellar streams in the Milky Way's halo, we hope 
that the result presented here will also be useful as a first step in the 
numerical studies of disrupting stellar substructures in the Galaxy. 

Acknowledgements 

We thank James Binney and Martin Smith for providing comments 
on a draft version of this paper. We also thank the referee, Paolo 
Miocchi, for helpful comments. SJ would also like to thank Eric 
Bell for discussions which inspired this work. 



REFERENCES 

Belokurov V., Zucker D. B., Evans N. W., Gilmore G., Vidrih S., 
Bramich D. M., Newberg H. J., Wyse R. F. G., Irwin M. J., Fell- 
hauer M., Hewett P. C, et al. 2006, ApJ, 642, L137 

Binney J., 2008, MNRAS, 386, L47 

Choi J.-H., Weinberg M. D., Katz N., 2007, MNRAS, 381, 987 
Dehnen W., OdenMrchen M., Grebel E. K., Rix H.-W, 2004, AJ, 

127, 2753 

Ibata R. A., Gilmore G., Irwin M. J., 1994, Nature, 370, 194 
Ibata R. A., Wyse R. F. G., Gilmore G., Irwin M. J., Suntzeff N. B., 

1997, AJ, 113, 634 
Jin S., Lynden-Bell D., 2007, MNRAS, 378, L64 
Jin S., Lynden-Bell D., 2008, MNRAS, 383, 1686 
Johnston K. V., Law D. R., Majewski S. R., 2005, ApJ, 619, 800 
Martin N. R, Jin S., 2009, ApJ, submitted 
Montuori M., Capuzzo-Dolcetta R., Di Matteo P., Lepinette A., 

Miocchi P, 2007, ApJ, 659, 1212 
Odenkirchen M., Grebel E. K., Dehnen W., Rix H.-W., Yanny B., 

Newberg H. J., Rockosi C. M., Martmez-Delgado D., Brinkmann 

J., Pier J. R., 2003, AJ, 126, 2385 
Odenkirchen M., Grebel E. K., Rockosi C. M., Dehnen W., Ibata 

R., Rix H.-W., Stolte A., Wolf C, et al. 2001, ApJ, 548, L165 
Paczynski B., 1990, ApJ, 348, 485 



